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We study bosons loaded in a one-dimensional optical lattice of two-fold p-orbital degeneracy at 
each site. Our numerical simulations find an anti-ferro-orbital Pcc+ipy, a homogeneous p x Mott in- 
sulator phase and two kinds of superfluid phases distinguished by the orbital order (anti-ferro-orbital 
and para-orbital). The anti-ferro-orbital order breaks time reversal symmetry. Experimentally ob- 
servable evidence is predicted for the phase transition between the two different superfluid phases. 
We also discover that the quantum noise measurement is able to provide a concrete evidence of time 
reversal symmetry breaking in the first Mott phase. 
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The last decade has witnessed rapid experimental de- 
velopments in the studies of ultracold quantum gases in 
optical lattices [J-Q. Ultracold bosons and fermions in 
optical lattices provide robust and controllable systems 
to study correlated quantum condensed matter physics 
beyond the scope of conventional many-body systems. 
For example, quantum phases and quantum phase tran- 
sitions of bosons in optical lattices were studied in recent 
years HQ. 

Besides exploring quantum phases of bosons in the low- 
est band, the possibilities of observing exotic phases of 
higher orbital band bosons were also put forward [H, ^ . 
In more recent experiments, a superfluid phase with a 
complex order in a two dimensional p-orbital band bo- 
son system fTol-[l2| has been reported, and evidences 
of exotic f-band superfluid phases have also been ob- 
served [13j- Theoretically exotic phases predicted for 
both fermions EMI and bosons [19-26] with orbital 
degrees of freedom are attracting considerable attention. 
A phase-sensitive scheme of detecting the complex or- 
der of the p^+ip^ superfluid (SF) was proposed very re- 
cently H3. 

In this article, we study bosons loaded in p x and p y 
orbits of a one dimensional (ID) optical lattice at zero 
temperature with both the numerical simulations and 
field-theoretical methods. We find two SF phases distin- 
guished by an orbital order — an anti-ferro-orbital (AFO) 
SF and a para-orbital (PO) SF, and two Mott insulating 
phases — an AFO Mott and a p x Mott phase (FIG. [1]). 
The AFO order is a staggered orbital current (p x ± ip y ) 
order [H, S3- In the AFO SF phase, the inter-band 
phase difference is locked at ±^ and the spontaneous 
AFO (p x ±ip y ) order in this phase breaks the time rever- 
sal symmetry (TRS), whereas the fluctuations of the rel- 
ative phase restore the TRS in the PO SF phase. Based 
on our results, we propose an experimental method to 
distinguish different phases by measuring momentum dis- 
tribution (FIG. [2]), instead of directly measuring the local 
current flow resulting from TRS breaking. In this way the 
PO to AFO quantum phase transition, associated with 



TRS breaking, can be observed in experiments. The fi- 
nite momentum peaks in the momentum distribution of 
the AFO SF phase make it distinguishable from the con- 
ventional ID SF phases HEHl]. In the AFO Mott phase 
the quantum noise measurement will be able to provide 
a concrete evidence of spontaneous TRS breaking. 

Experimental proposal. — The system we shall propose 
is a ID lattice elongated along the x direction, and each 
lattice site has a rotation symmetry in the x-y plane 
(FIG.[T]). In other words, the p x andp y orbits are locally 
degenerate, but the hopping differs significantly. Such 
a ID system can be realized from a 2D optical lattice. 
Suppose the 2D optical lattice is formed by different laser 
beams in the x and y directions and the lattice poten- 
tial reads V = V x ^>\v?{k x x) + V y sin 2 (k y y), where V x and 
k x (V y and k y ) are the strength and wave numbers of 
the laser beams in the x (y) direction. In tight binding 
approximation, we can use the harmonic wavefunctions 
to approximate the Wannier functions. In the harmonic 
approximation, the local isotropy (rotation symmetry of 
each site in the x-y plane) requires V x k 2 = V y ky. This re- 
lation, which guarantees the (approximate) two-fold or- 
bital degeneracy at each lattice site, can be somewhat 
surprisingly well held in the ID limit by taking the lattice 
potential depth V y ^> V x and simultaneously the lattice 
constants a y (= ■^ L ) > a x (= ^ L ). As a result, the local 
isotropy is maintained, but the system has stronger po- 
tential and larger lattice spacing in the y direction than x, 
which makes the hopping in the y direction much smaller 
than that in the x direction. As an example, we take 
V x /E RjX = 6, V y /E Rj y = 24 and a y /a x = V% where 
Er^q, = h 2 k^/2m is the recoil energy in the a direction. 
Such a system is locally isotropic, while the hopping of 
p y orbital in the y direction is smaller than one percent of 
the hopping of p x orbital in the x direction, and the sys - 



tem is dynamically ID. With the technique in Ref. 



bosons can be loaded in the p-orbits of optical lattices. 

The life time and the phase coherence of p-orbital 
bosons could be enhanced by applying double wells as 
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FIG. 1. (Color online). Phase diagram of a one dimensional 
lattice Bose gas with p x and p y orbital degrees of freedom. 
The upper panel shows the sketch of experimental setup we 
proposed. The green circles are used to denote the require- 
ment of the approximate local isotropy of the lattice potential 
at each site. The lowest Mott lobe (with filling v — 1) is 
dominated by p x bosons. The Mott state (with v > 1) has an 
AFO order (see text). We do not claim another phase for the 
tiny tip of the second Mott lobe beyond the red line because 
of numerical errors. For sufficiently large hopping t x or for 
low filling, the Bose gas has a crossover from PO SF to a p x 
SF phase, which will not be discussed in this work. 



in the experiment . The double well lattice gives un- 
equal band gaps (between the p and the lower and higher 
orbital bands), suppressing the decay of p-band bosons. 
Here, for the sake of illustrating the salient features of 
the ID degenerate p-orbits, we employ a simple single- 
well lattice to simplify analysis. Given the consistency of 
the experiment with the theory on a simple square 



lattice Jjl, |20] , we do not expect the complexity induced 
by double wells would alter the basic properties and the 
understanding of the reported quantum phases. 

In the ID limit, the Hamiltonian describing bosons 
loaded in these p x and p y orbits reads [20| 



<jj'> 

Z^j 2 



-t x al(j)a x (f) - 
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Here a x (j) (a y (j)) is the annihilation operator for p x 
(p y ) orbital at site j. The discrete variable j labels 
the sites of the ID chain, with the lattice constant a x . 
The local particle number operator h(j) is denned as 
^2 a=x y ^a(j)^a(j): and the local angular momentum op- 
erator L z (j) is defined as o e a ^(— ia^ where 
the superscripts a and j3 run over x and y. U (> 0) is the 
repulsive Hubbard interaction. The average number of 
bosons per site is fixed by chemical potential \i. t x (< 0) 
is the longitudinal hopping of p x bosons, and t y (> 0) 
is the transverse hopping of p y bosons (FIG. [1]). Due to 



FIG. 2. (Color online). Momentum distributions n(k). (a) 
shows the ID momentum distribution hid(k x ) in a geometri- 
cally 2D experimental setting (see paragraph " Experimental 
proposal") for PO SF phase. p x peaks at k x a x = (2j + 1)tt (j 
is some integer) are sharp, and p y peaks at k x a x = 2 jit are 
broad. The insets show that the double logarithmic plot of 
nid(k x ) near the sharp (broad) peaks is linear (non- linear), 
(b) shows the sketch of 2D momentum distributions (n(k)) in 
different phases (PO SF, AFO SF and AFO Mott from right 
to left). In three subgraphs the horizontal (vertical) axis is 
k y a y /7T {k x a x /n). The purple wiggles along each subgraph 
shows hid(k x ). In the AFO SF phase, the p y peaks which are 
broad in PO SF, are replaced by sharp peaks. In the AFO 
Mott phase, there are no sharp peaks. 



anisotropy of the p-orbits, the longitudinal hopping ("a 
bond" ) is typically much larger than the transverse hop- 
ping ( u 7r bond") [19l42l|. In this paper the ratio \t x /t y \ 
is taken as 9, which corresponds to V x « 6Er iX (based 

on the estimation \t x /t y \ « ^i^ 2 \J~^~ ~ 6) under tight 

binding approximation) . 



The hopping term has a U(l)xU(l) symmetry, which 
is a a (j) [e™° u e™*%pap(j) : with a 1 ° 
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It appears that the particle num- 



bers of p x and p y components, N x — ^2jn x (j) and 
N y = Y^j n y{j) are separately conserved. However the 
pair hopping term a y a y a x a x from L 2 Z does not conserve 
N x and N y separately, and thus breaks the U(l)xU(l) 
symmetry. Only the total particle number N = N x + N y 
is conserved. The U(l)xU(l) symmetry is reduced to 
U(l)xZ 2 denned as a a (j) -> [e ia ° e e i(7z %] a pap(j). 

Numerical Methods. — We use a matrix product state 
to represent the ground state. The ground state is ob- 
tained by iterative optimization [30|. White's correction 
is implemented to avoid potential trapping in the iter- 
ation procedure [31]. An open boundary condition is 
adopted in this work. The good quantum number we 
used in our numerical calculation is the total particle 
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number N. The chemical potential is calculated as the 
energy it takes to add a particle (hole) to the many-body 



state 



32, 3 



The Mott gaps are determined by extrap- 
olating to the thermodynamic limit. The largest sys- 
tem studied has 120 sites, which is large enough to com- 
pare with the experiments on ID quantum gases. With 
the numerical method the ground state phase diagram 
of Hamiltonian in Eq. (pQ) is mapped out and shown in 
FIG.[U The phase boundary of the Mott insulating phase 
is determined by the vanishing of Mott gap. The phase 
boundary between the AFO and PO SF phases is de- 
termined by the vanishing of the Z2 order parameter, 
defined as L z = J^ . (e l ® J L z (j)), with Q = tt. We use a 
system of 40 sites to determine this phase boundary. The 
existence of the AFO and PO SF phases is verified for a 
system with up to 100 sites. The central focus of this pa- 
per is the finding of unexpected quantum orbital phases 
in a ID optical lattice. A more accurate calculation of 
phase boundaries is left for future study. 

Mott phases.— For the Mott phases (FIG. [Q), the 
filling factor v — (n(j)) at each site is commensurate. 
The occupation number for each orbit, both (a x a x ) and 
(a y a y ) are incommensurate. For filling v greater than 1, 
the Mott phase features a complex order (a x (j)a y (j)) ~ 
e *Qi e <f with f = zb spontaneously chosen, which breaks 
the U(l)xZ2 symmetry down to U(l). Equivalently this 
Mott state has a staggered angular momentum order 
(L z (j)) ~ e 1 ® 3 ' . The order parameter L z is finite. With- 
out loss of generality, we have assumed L z is positive. 
This Z2 order also breaks the TRS, because finite L z 
means a finite local vortex-like current flow. For filling 
v equal to 1, the Z2 order does not exist for \t y \ <C \t x \. 
We call it p x Mott since p x boson dominates this Mott 
phase, i.e., (a x a x ) > (a y a y ) for \t x \ > \t y \. 

AFO superfluid phase. — By increasing the hopping 
the system goes into the SF phase when the Mott gap 
closes. The system has a phase transition from the AFO 
Mott phase to the AFO SF phase (FIG. [Q). Since the 
Z2 symmetry is broken in this SF phase, it behaves like 
a single component SF phase far from the Z2 critical 
point. This AFO SF phase is thus characterized by an 
algebraic correlation (aj(j , )af (j)) ~ \j — j 

a\(j) = e^at(j) + i^oj, (j), with £ € (0, 1]. £ 
limit of t x /U 0. The phase correlations of original bo- 



K l 2 , where 
1 in the 



son operators, defined as G a j3(j,j f ) - 
given by G xx (j,j') ~ e^-^\j - f 

ip iQj\n-i'\~K/2 and q 



{<{j)ap{j')), 



are 



' K/2 , G x (j,f) ~ 

ie^\j- j'\~ K/2 and G yy (j,f) ~ \j-j'\- K ' /2 in the AFO 
SF phase. We emphasize that the TRS is broken in this 
phase, because the off-diagonal correlation G xy (j,j f ) is 
complex. The key feature is that the power law decay 



-K/2 



) correlations (G xx G xy and G yy ) exhibit 



(li - 

the same power exponent K/2. In this phase, the rela- 
tive phase between the p x and p y SF components is 
locked. The two components share the same U(l) phase 
(p+ at low energy. The Lagrangian describing phase flue- 
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FIG. 3. (Color online). Properties of the AFO to PO phase 
transitions with total filling v — 1.5. (a) shows the Z2 order 
parameter L z . Our results indicate a continuous phase tran- 
sition of the AFO order, (b) shows the filling of p x and p y 
bosons. The p y component does not vanish across the phase 
transition, (c) shows the ratio K x /K y in the PO SF phase. 



tuations is 



1 

2ttK 



[v+\d T (p+) 2 + v + (d x (p + ) 2 



(2) 



The Bose liquid is completely characterized by the sound 
velocity v+ and Luttinger parameter K. For K < ^, the 
AFO SF phase is stable against the periodic lattice po- 
tential; for K > ^ at commensurate filling, this phase is 
unstable and undergoes a localization transition towards 
the Mott phase [11 EH- The AFO order is preserved 
across the localization transition in our system. 

PO superfluid phase. — The AFO order disappears for 
larger hopping, and the AFO SF gives way to the PO 
SF. The behavior of the Z2 order parameter L z and oc- 
cupation numbers of p x and p y bosons across the phase 
transition is shown in FIG. 03 The Z2 order is de- 
stroyed by quantum fluctuations of cp- and thus the 
TRS is restored in the PO SF phase. The phase cor- 
relations of the original bosons in this phase are given 
by G xx (j,j>) ~ e*OW-i')|j _ j>\-K./2 an d G yy (j,j') ~ 
\j ~ j'\~ K y/ 2 . The phase coherence between p x and p y 
components — G xy (j,j') — vanishes in this phase. By nu- 
merical simulations, we find K x <C K y . Following Hal- 
dane [36], an analytic expression estimating K x /K y is 

derived, ^ ~ ^/||^-|, where n x (n y ) is the filling of 

Px (Py) bosons. This is in qualitative agreement with the 
numerical results (FIG.[3fc)). The Lagrangian describing 
this PO SF phase is 



+\(d T (p x )(d T <p y ), 



[v^idrifa) 2 + V a (d x (f a ) 2 



(3) 



where cp x ((p y ) is the phase of p x (p y ) SF component. The 
mixing term (A) is much smaller than the kinetic term 

(2^) in our s y stem - 

Quantum phase transition from AFO to PO in the su- 
perfluid phases. — The phase transition from the AFO 
SF to the PO SF is described by a sine-Gordon model of 
the relative phase cp- . The Lagrangian is 



C[<p. 



2irK- 



[v^tdrip-) 2 + v_(^(^_) 2 ]+mcos(2(^_), 

(4) 
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where rn is estimated as m « \Un x n y . When rn is 
greater than some critical m r (K_), the sine-Gordon the- 
ory is in a gapped phase [37], and therefore ip _ field is 
locked at one minimum of racos(2<^_). Such an orbital 
gapped phase is the AFO SF phase. When m < m c (K-), 
the sine-Gordon term mcos(2(/?_) is irrelevant in the 
sense of renormalization group, and the theory is in a 
gapless phase for which ip_ is unlocked. This orbital 
gapless phase is the PO SF phase. We emphasize here 
the sine-Gordon term is not perturbative in our model. 

Experimental signatures. — Since the quantum phases 
we have found are characterized by distinct phase cor- 
relations, measuring the momentum distributions by 
time-of-flight (TOF) will distinguish different phases. 
Assuming the interaction of atoms in TOF is negli- 
gible, the density obtained in TOF experiment (ntof) 
measures the momentum distribution given as n(k) oc 
l^a=x,y^nGaa(k x a x + 2n7r)</>* (k)0 a (k) in our system 
(shown in FIG. [2]), with G aa (k x a x ) the Fourier transform 
of G aa (j,j') and a (k) the momentum-space form of p- 
orbital Wannier functions. The ID momentum distribu- 
tion is defined as n\d{k x ) = J dk y h(k). The strong peaks 
of momentum distribution at finite momenta k x = tt / 'a x 
mod 2ir/a x will distinguish the AFO SF phases from con- 
ventional ID SF phases j^ESt The Luttin ger parame- 
ter can be estimated by measuring momentum distribu- 
tion since it behaves as \og(hid(k x )) ~ (4f — 1) \og(\k x ± 
7r/a x \) for k x near peaks ±7r/a x . Including the effect 
of Wannier functions ^(k), which are smooth slowly 
varying functions, does not alter appreciably the form of 
nid(k x ) near ±it/a x . The effects of the harmonic trap- 
ping potential [38| and the interactions in TOF, which 
cause difficulties of extracting the Luttinger parameter, 
are left for future studies. 

The AFO order in the Mott phase will have exper- 
imental signatures in the quantum noise measurement. 
The quantum noise is defined as C(d) = J <i 2 Rg(R, d), 
with flf(R, d) = (n to f (R + ~d)n to f(R - |d)) - (n to f(R + 
|d)) (ntof (R-|d)), and R= (R x ,R y ),d= (d X: d y ). The 
brackets ( ) denote statistical averages of independently 
acquired TOF images in experiments [39|. For the Mott 
phases in our proposed 2D optical lattice, g(R, d) is given 

by 

<?(R, d) = L {Ek {ft d ~ K ) (Cxxn x + (yyUyf 

+ Ek^ (2) &-K + Q x + Q y ) 

x \CxyGxy + CyxG X y\ } j (5) 

where t is the time of flight and L is the number of 
lattice sites. Here, Q xy = G^(0,0), ( a /3 ~ (R a + 
\d a ){R - \d p ), Q, = (£,0), Q y = (0,£) and 
K = 2jiQ x + 2j 2 Q y (ji and j 2 are integers). In Eq. (j5j) 
the smooth Gaussian part of Wannier functions ^^(k) 
is approximated by a constant function, which is typical 
in quantum noise measurement [39|. The center of the 



trapped cold gas is taken as the origin of coordinates here. 
The sharp peaks of C(d) at d = d = ^ (K - Q x - Q y ) 
signify that the off-diagonal term Q xy is finite, which dis- 
tinguishes the AFO Mott state from the p x Mott. The 
experimental signature of an imaginary Q xy is predicted 
to be that g(R, d) exhibits nodal lines at R || do. Q xy be- 
ing imaginary tells a local vortex-like current flow, which 
is a concrete evidence for the TRS breaking. 
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